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1. Introduction 

Integrable operators with kernels of the form 

W(x,y)= nx)9{y) - ny)9ix) (1.1) 
x-y 

have applications in quantum field theory and random matrix theory, where they are used 
to describe the asymptotic distribution of large random matrices; see [5, 17, 18]. Tracy 
and Widom [19] observed that many important distributions in random matrix theory can 
be defined using solutions of systems 
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where m(x), a(x), (3(x) and 7(2;) are polynomials. In [17, 18], Tracy and Widom considered 
the Airy and Bessel kernels which describe the soft and hard edges of generalized unitary 
ensembles, and proved the apparently miraculous identities that the operators with these 
kernels were squares of self-adjoint Hankel operators; then they used this property to 
compute their eigenf unctions and eigenvalues. 

Let if be a separable Hilbert space, and L 2 ((0, 00); dx; K) be the Bochner-Lebesgue 
space of strongly measurable functions (ft '■ (0, 00) — ► K such that f Q || (j)(x)\\ 2 K dx < 00. 
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In [2] we considered a general class of differential equations which gives rise to integrable 
operators that have the form W = T*T : where V : L 2 (0, oo) — > £ 2 ((0, oo);K) is a continu- 
ous Hankel operator. The general theorem of [2] specialised to the Airy and Bessel kernels, 
and in this paper we prove related results which deal with other integrable operators. 

After recalling some definitions, we state and prove the main Theorem 1.2, and then 
in section 2 give some applications. In section 3 we discuss the scope of Theorem 1.2 as 
it applies to (1.2) in reduced forms, and then in section 4 prove another theorem which 
encompasses other applications, as in section 5. 

Definition (Integrable operators) Let I be a subinterval of R. An integrable operator on 
L 2 {I;dx) is a continuous linear operator W with kernel 

W(x,y) = 2j2 fjiX)gjiy) (z,yel;z*y) (1.3) 
A — ' x — y 

where fj, gj are bounded and measurable functions such that Y^=i fj( x )9j( x ) = almost 
everywhere on I. 

Lemma 1.1. Suppose further that the fj and gj are real-valued. Then W is self-adjoint 
if and only if 

w{ ) = {M^m (XtyeI , x ^ v) (L4) 

x-y 

where v(x) = col[/i(x), . . . , /„(x); gi(x), . . . , g n (x)] and 



J = 



-In 
In 



(1.5) 



with identity matrix I n € M n (R) and the usual inner product on R 2n . 

Proof. Clearly the kernel of W is symmetric if and only if the numerator of W(x,y) is 
skew-symmetric, in which case we can write 

n n n 

2Y,f^)9M = T,M x )9j(y)-T,fM9j^)^ (1-6) 

i=i j=i j=i 

and the matrix expression follows directly. 

□ 

Definition (Hankel operator) Let <fi G L 2 ((0, oo); K). The Hankel operator with symbol 
(f) is the integral operator 

/■OO 

T 4> f{s)= / <P(s + t)f(t)dt. (1.7) 
2 



Nehari's theorem [13] gives a sufficient condition for : L 2 (0, oo) — > L 2 ((0, oo); K) to be 
continuous and gives an expression for the operator norm ||r||; while the Hilbert-Schmidt 
norm of satisfies 

/•OO 

\\^4hs= / sU(s)\\ 2 K ds (1.8) 



JO 

when y/s(j)(s) G £ 2 ((0, oo); K). Clearly is self-adjoint when K = R. 

Given an integrable operator, it is often valuable to identify a Hankel operator 
such that T^r^ = W and to determine whether W is of trace class. In particular, when 
K = R and W = the spectral resolution of the self-adjoint operator T^, determines 
the spectral resolution of W. This is the basis of the successful calculations in [17, 18, 
19], which also exploited the fact the eigenvectors of can be relatively easy to analyze. 
In [12], Megretskh, Peller and Treil characterized the spectral multiplicity function of a 
self-adjoint Hankel operator. Further, in applications to determinantal point fields as in 
[15], one often wishes to show that W is of trace class and satisfies < W < I. 

We consider first the case associated with the differential equation 

du 

— = J(Oix + O + tt-ix~ 1 )v. (1.9) 
ax 

Theorem 1.2. Suppose that Qi,Qq and 0_i are reai symmetric (2n) x (2n) constant 
matrices such that fii > and — > 0. Suppose further that v satisfies (1.9), and that 
v(x) and v(x)/x are bounded functions in L 2 ((0, oo); dx; R 4n ). Then there exists a real 
linear subspace K of R 4n with dim(K) < rank(O x ) +rank(0_i) and (p £ £ 2 ((0, oo); dx; K) 
such that : L 2 (0, oo) — > L 2 ((0, oo); dx; K) is continuous and the kernel 

W{x,y)= {Mx)Mv)) (1.10) 
x-y 

factors as 

w = r;r 4> . (1.11) 

In particular, if rank(Oi) + rank(0_i) = 1, then W = T 2 for some i> G L 2 (0, oo). 

Proof. By Lemma 2.1 of [2], which essentially depends upon the continuity of the Hilbert 
transform on L 2 (R), we know that W is also a continuous linear operator on L 2 (0, oo). 
We have 
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where dv/dx and dv/dy satisfy (1.9). Now J 2 = —I and J* = — J, so the matrices involved 
in the differential equation satisfy 

j 2 (n lX + o + ft-ix' 1 ) + (o l2 , + o + o.iy- 1 ) r j 

= -Oix + Oiy-0_ix _1 + 0_iy _1 ; (1.13) 
so by dividing by x — y, we obtain 

(I + h) w(x ' v) = -<""'<*>■«'<''» + ("-^ • T 1 )' < L14 > 

We introduce the positive roots of the positive semidefinite matrices Oi and 0_i, and the 
column vector (f)(x) — col[y/Tliv(x), v 7 — 0,_ 1 v(x)/x]; so that, G £ 2 ((0, oo); da;; R 4n ) and 

idi + dy^) Wix ' y) = -ttWMv)) (x,y> 0). (1.15) 
Integrating this equation, we obtain 

poo 

W{x, y) = / (<j){x + t),<f>{y + t)) dt + g(x - y) (1.16) 
Jo 

for some differentiable function g; but W(x,y) and the integral converge to as x — > oo 
or y — > oo; so (7 = 0. Hence I4 7 = T^r^, and we deduce that is a continuous Hankel 
operator from L 2 (0, 00) -> L 2 ((0, 00); da;; R 4n ). 

Finally, we observe that 4> takes values in a linear subspace K of 
range(Oi) ©range(0_i) which has dimension less than or equal to rank(Oi) +rank(0_i). 

If the sum of the ranks of Oi and 0_i equals one, then takes values in a one- 
dimensional real linear subspace of R 4n , so <f>(x) = eip(x) for some unit vector e G R 4n 
and some ip G L 2 (0, 00); hence T^, is self-adjoint and W = T 2 . 

□ 

Remarks 1.3. (i) In their applications to random matrix theory, Tracy and Widom con- 
sidered integrable operators on unions of intervals such as U^L-Jg^-i, 0,2 j] - Many analytical 
problems reduce to considering one interval at a time, and so are addressed by the current 
paper. In a subsequent article [20], they generalized their results to kernels of the form 

(C<p(x),<p(y)) 
x-y 

where C is an antisymmetric n x n matrix. 

(ii) Theorem 1.2 has an analogue for discrete kernels on Z + , as in [11]. 
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2. Factorization for some differential equations with simple poles 

In this section we consider how Theorem 1.2 applies to some differential equations that 
are satisfied by familiar special functions. 

2.1. The Airy equation 

The Airy function Ai gives rise to a solution u(x) = Ai(x + s) of the differential 
equation 

u" (x) = (s + x)u(x) (x > 0). (2.1) 
Since the standard asymptotic formula for the Airy function [16, p. 18] gives 

Ai(x) = ^T/^(l + 0{x~^) exp(-^ 3 / 2 ) (x - oo) (2.2) 

the hypotheses of Theorem 1.2 are satisfied by v(x) = col[<j)(x), <f>'(x)] where <f>(s) 

= Ai(x + s); so W = where is a Hankel operator in the Hilbert-Schmidt class. See 

[2, 17] for more details. 

2.2. The Laguerre equation 

The Laguerre equation [14] may be expressed as 



u"(x)+(-\ + ^±±)u(x) = 0, 



(2.3) 



with solution u(x) = xe x / 2 Ln\x) where 



L *i x) = x -^^ xn+le "1 {x>0) 



(2.4) 



is the Laguerre polynomial of degree n and parameter a = 1. The Laplace transform of u 
is the rational function 



C(u; A) = (n+1) 



> -1/2). 



Theorem 1.2 applies directly to the system 



d 


u(x) 




dx 


u'{x)_ 





1 

l/4-(n+l)/x 



u(x) 
u'{x) 



and gives the formula 

u'(x)u(y) - u'(x)u(y) 
x-y 



= (n+l) 
5 



1 u(x + t)u(y + t) 
o (x + t)(y + t) 



dt 



(2.5) 



(2.6) 



where <f>(x) = u(x)/x gives a Hankel operator of Hilbert-Schmidt type. 

2.3. The Bessel equation 

The differential equation 

1 





u"(x) H — 
x 



(2.7) 



has solution it (a;) = y/xJi(2y/x), where J\ is the Bessel function of the first kind of order 
one. The Laplace transform of it satisfies C(u, A) = A -2 exp(— 1/A). Now the standard 
asymptotic formula for the Bessel function [10, p. 171] shows that 



Mx) = x -^x~ 3/4 cos [2y/x- — (x^oo); 
a; \/7r V 4 / 



(2.8) 



so (p belongs to L 2 (0, oo); hence one can follow the proof of Theorem 1.2 and derive the 
formula 



u(x)u'(y) — u'(x)u(y) 



' u(x + t)u(y + t) 
o ( x + t)(y + t) 



dt 



x-y 

Here the Hankel operator is not Hilbert-Schmidt. 

2.4. The Carleman operator with multiple spectrum 

The system 



(x,y > 0). 



(2.9) 



d 


7" 




"0 1/x' 




7" 


dx 
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(2.10) 



has the form considered in Theorem 1.2 and evidently has solution f(x) = logx and 
g(x) = 1; further 



W(x,y) = 



log x — log y 

x-y 



dt 



o (x + t)(y + t) 



(x,y> 0) 



(2.11) 



has a similar form to an integrable kernel, except that log a: is unbounded. Power showed 
that Carleman's operator T, where 



Fh(x) = [ 
Jo 



h(y)dy 

x + y 



(h E L 2 (0,oo)), 



(2.12) 



is continuous on L 2 (0, oo) and has spectrum [0, 7r] with spectral multiplicity two; see [13]. 
Hence W is a continuous linear operator on L 2 (0, oo) with spectrum [0, 7r 2 ] with multiplicity 
two. This example illustrates that simple differential equations can give positive definite 
Hankel operators with multiple spectra. 
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and let 



2.5. Parabolic cylinder functions: non factorization 

For $tp > — 1, let D p be the parabolic cylinder function, which satisfies 

DJx)D'(y) - D'(x)DJy) 

H(x, y) = PK 1 pK } pK J p . (2.13) 

x-y 

Then ±H is not the square of a self-adjoint Hankel operator T^. By following the proof 
of Theorem 1.2, we obtain 

(a± + Y y ) H{x ' y) = -2 {x + ( 2 - 14 ) 

where (x + y)D p (x)D p (y)/2 cannot equal ±4>(x)(/)(y)] indeed, for suitable xi,X2 > 0, the 
2x2 matrix [xj + Xk]j,k=i,2 has both positive and negative eigenvalues. When n is a 
nonnegative integer, D n is known as a Hermite function, and may be written 

Ux) = (n!)- 1 /2( 27r) -V4 ( _ 1) n e xV4A_ e - 2 /2. (2 . 15) 

Aubrun [1] considers self-adjoint Hankel operators T ( j }n and T ( j )n+1 such that 
H = (l/2)(T ( f >n r < f >n+1 + r^^r^J; this gives information about the singular numbers of 
H. For Hankel squares one has more precise information about the eigenvalues. Borodin 
and Okounkov [4] have considered the discrete Hermite kernel, and derived the formula 



<fim+l(s)<fin(s) ~ (pm(s)(pn+l(s) 



m — n 



poo 

/ 4> m (t)4> n (t) dt (m,?j = I l,... I ffl/n). (2.16) 

J s 



Here the variable n is the degree of the Hermite polynomial factor in <fi n , and (2.16) is 
essentially different from (1.11). 

3. Reducing to standard form 

Definition {Operator monotone) Let / be an interval in R. A continuous function to : I — > 
R is operator monotone increasing if, whenever S and T are continuous and self-adjoint 
linear operators on Hilbert space that have spectra in /, 

S <T ^u(S) <u(T). (3.1) 



Further, u> is operator monotone if and only if the matrices 

-uj(xj) - iv(x k )- 



Xj-Xk J j,k=l,...,m 
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(3.2) 



with diagonal entries u'(xj), are positive semidefinite for all m = 2, 3, . . . and Xj G I with 
j = 1, . . . , m. In [2] we used Loewner's characterization of operator monotone functions, 
which shows in particular that an operator monotone function on (0, oo) extends to an 
analytic function on a domain U containing (0, oo) as in [9, p. 541]. 

Theorem 1.2 shows that, under mild technical conditions, W admits of a factorization 
W = r^r^ whenever uj(x) = (0(a;)£,£) is operator monotone on (0, oo), for all £ G R 4n 
where Q(x) = flix + Oo + fi-ia; -1 . 

Suppose that / has as an endpoint, and let U be a domain that contains I. Suppose 
that A(z) is a matrix function into M2 n (C) that is analytic on U, except for an isolated 
singularity at z = 0, and that 

flu 

^=A(z)v (zeU). (3.3) 
dz 

By a standard change of variable, we mean w(z) = T(z)v(z), where the analytic function 
T : U — > M n (C) has T(z) invertible as a matrix for each z G U. The following result is 
commonly known as Birkhoff 's normal form, although the first correct statement and proof 
is due to Turrettin [21]. Gantmacher considered some related examples which resemble 
(2.10) in [6, p. 146]. 

Proposition 3.1. Suppose that A{z) = J2k=-oo^kZ k is a Laurent expansion that con- 
verges for all z 7^ 0. Then there exists a standard change of variable that reduces (3.3) 
to 

rfj = hr + — r- (3 - 4) 

Further, if the eigenvalues Xj have differences Xj — Xk that are never equal to a natural 
integer, then one can remove the term in A_2- 

The appearance of the term A_ 2 in (3.4) is important, since — 1/x 2 is not operator 
monotone on (0, oo) by [9, p. 554]. So we cannot simply adapt the proofs in section 1 to 
deal with the case in which A_ 2 appears. However, if u> is operator monotone on (0, oo), 
then uj(^/i) is likewise. This suggests the change of independent variable x = Vt, which 
we exploit in the examples in section 5. Further, we adjust the definition of the kernel and 
the Hankel operators so that we can obtain a factorization theorem in the next section. 

Definition (Hankel operator) For I = (1, oo), we use the Hankel operator 

l>0(a;) = / ij(xy)g(y)— (3.5) 

j\ y 



where ip G L 2 ((l, oo); dy/y; K); whereas for / = (0, 1) we use 

I>0r)= f p(xy)h(y)^ (3.6) 

Jo y 

where p G L 2 ((0, 1); dy/y; K). These definitions reduce to the case / = (0, oo) in section 1 
by the changes of variables y = e l and y = e~ l respectively. 

4. Factorization theorem for differential equations with double poles 

In this section we consider the differential equation 

x— = JiQxx + Qo + aJ + Q- 1 x~ 1 )v. (4.1) 
dx 

To accommodate forthcoming examples, we have introduced the skew-symmetric matrix 
a J into the constant term for some a G R. 

Theorem 4.1. Suppose that Oi,O and are real symmetric (2n) x (2n) constant 
matrices such that fii > and — > 0. Suppose further that v satisfies (4.1), and 
that x a v(x) and x a ~ 1 v(x) are bounded functions in L 2 ((l, oo); dx; R 4n ). Then there ex- 
ists a real linear subspace K of R 4n with dim(K) < rank(Oi) + rank(0_i) and <j> G 
L 2 ((l, oo);dx/x; K) such that : L 2 ((l, oo);dx/x) — * L 2 ((l, oo);dx/x; K) is continuous 
and the kernel 

W(x, y) = ™— (Jv(x),v(y)) (4.2) 

x y 

factors as 

w = r;rv (4.3) 

In particular, if rank(Oi) + rank(0_i) = 1, then W = T 2 b for some ip G L 2 ((l, oo); dx/x). 
Proof. We observe that by homogeneity 

/ d d\{xyY 2a+l ^ 2 n {xyY 2a+l ^ 2 , n , . 

( x ir + yir) = 2a^ (x,y>0;x^y), 4.4 

V ox oyJ x — y x — y 

and hence 

9 , d \„ r , , n (xy)&+W 



/ O O \ [XV) y J/ 

\ X di +V di) W(x - v) = 2 " x-v ( JV(XI 



+ <s^ ((j *^ )> + (j ^ * >}i (46) 



a; — j/ V\ da; ' / \ ' dyl 
where the matrices involved in the final terms in (4.5) are 

J 2 {$l x x + Q + aJ + 0_ia; _1 ) + (ttiy + tt - aJ + 0_iy _1 ) J* J 
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= -Qx(x -y) + Cl-i(x - y)/{xy) - 2aJ. 
By cancelling the terms that involve J, we obtain 
d d 



(x—+y—)w(x, y) = -(xy) 2a+1 ^ 2 (n 1 v(x),v(y)} + (xy^+^^v^ v(y)). (4.6) 



dy 

We introduce the column vector 



y/n~ lX ^ a+1 ^ 2 v{x) 

A /=n_ lX (2a-l)/2 u ( a .) 



(4.7) 



ip(x) ■ 

which belongs to £ 2 ((1, oo); dx/x; R 4n ) and satisfies 

( X lL + V ^y) Si ^ ty ^ T = ^(^))- 

Hence 

W(x,y) = j (iP(tx),i;(ty))j + h(x/y) (4.8) 

where h(x/y) — > as x — > oo or j/ — > oo; so /i = 0. One can conclude the proof by arguing 
as in Theorem 1.2. 

□ 

We consider later some examples in which 0_i =0. In this case, we can invoke the 
following existence theorem for solutions. 

Proposition 4.2. Suppose that the residue matrix A_i has eigenvalues Xj such that the 
differences Xj — Xk are never equal to a natural integer. Then the differential equation 

z^-X = (A z + A- 1 )X (4.9) 

ciz 

with X(z) G M 2n (C) has a non-trivial solution of the form X(z) = Y(z)z A - 1 , where Y is 
an entire matrix function of order one. 

For a proof see [8] , where Hille also discusses the asymptotic form of the solutions in 
terms of the Laplace transform. Note that the Laplace transform of (4.9) has a similar 
form to (4.9) itself and in particular has the residue matrix A_\ + I. 

5. Examples of factorization for differential equations with double poles 

5.1. Modified Bessel functions 

For < v < 1, MacDonald's function is defined by 

/>oo 

K v (z)= e- zcosht cosh(ut)dt (3fc>0) (5.1) 
Jo 
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and satisfies the modified Bessel equation z 2 K" + zK' v — {y 2 + z 2 )K v = 0; hence u{x) = 
y/xK v {2y/x) satisfies 



,1 v*-\ 

U (x) = [ 1 — rr- 

\x Ax 2 



u(x). 



(5.2) 



By [7, 8.451], K v (x) decays exponentially as x — > oo. We can apply Theorem 4.1 to the 
system 

d r,, i r _i 1 1 r„ i 

(5.3) 



d 


u 




-1 1" 




u 


dx 


w 




x-2 + \{v 2 -1) 2_ 




w 



so that, in terms of Theorem 4.1, 



1 




O = 



-2 + l(^ 2 -l) 

3 
2 



-1 



, a = -1/2, O_ 1 = 0, (5.4) 



where the residue matrix A-i = JOo — al has eigenvalues (1/2) ± (v/2). 
Thus one obtains 



u(x)v(y) — u(y)v(x) 
x-y 

and after some reduction, one deduces that 

K v {2^)^yK' v {2^y) - ^K' v (2^)K v (2Jy) 



u{tx)u{ty) 



dt 



/oo 
K u {2yftx)K v {2s/ty)dt (5.5) 



x-y 

where the right-hand side is the square of a Hankel operator of Hilbert-Schmidt class. 

5.2. Bessel functions 

The Bessel function J v satisfies x 2 J" + xJ' u + (x 2 — u 2 )J v = 0, and hence u = 
^fxJ v {2^fx) satisfies 



u"(x)+^+ 1 -^f)u(x) = 0. 



~x ' Ax 2 

One can apply Theorem 4.1 with some obvious sign changes to the system 



d 


u 




-1 


r 




u 


dx 


w 




-x-2-\{l-v 2 ) 


2 




w 



(5.6) 



(5.7) 



and after some reduction one obtains an identity from [18] 
xJl(2^)M2^y) - Jvpyffiy/yJlpy/y) _ " j 



x-y 



J v {2y/tx)J v {2y/ty)dt. (5.8) 
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5.3. Whittaker's functions 

The homogeneous confluent hypergeometric equation may be reduced to Whittaker's equa- 
tion 

and the solutions of this are known as Whittaker's functions. In [7, 9.227], the authors 
give a solution w(z) = W K ^(z) such that w(z) x e~ z l 2 z K as z — > oo along (0, oo). Kernels 
involving Whittaker's function appear in [3]. 

Proposition 5.3. Suppose that u(x) = W K>u (2y/x) for some k < 0. Then there exists a 
function $ e £ 2 ((1, oo);dx/x; K) for a separable Hilbert space K such that 



W(x,y) = (xy) 1 



/4 u{x)u'(y) - u'(x)u(y) 
x-y 



(5.10) 



factors as W = r|IV 

Proof. We can write, after a little reduction 



d 


u 


1 


-4 4" 




u 


dx 


V 


4 


x-2^v^- (i -z/ 2 ) -6 6 




V 



(5.11) 



which gives 



fl(x) 



1 


' x - 2ny/x- {\ - u 2 ) - 6 5 " 


1 


"0 -1" 


4 


5 -4 


~ 4 


1 



(5.12) 



hence a = —1/4. The function yfx is operator monotone increasing on (0, oo); indeed, 
4>t(x) = t 1 ' 4 /{t + x) belongs to L 2 ((0, oo); dt) with ||</>t||| 2 = n/(2y/x) and satisfies 



Vtdt 



x-y tt J (x + t)(y + t) 

1 f°° 
= - / (j) t (x)(j)t(y) dt. 
t Jo 

We observe that x 1 ^ 4: u(x)<pt(x) belongs to L 2 ((l, oo);dx/x) and since 

\dx + y d~) w ^ x,y " > = -i( x v) 1/4:u ( x ) u (y) 



(5.13) 



K 



+■ ^(xy) 1/4 u(x)u(y) J 4> t {x)4> t {y) dt 



(5.14) 
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we have 

W(x,y) = j f u(xs)u(ys)(xy) 1/4 -^ 



4J 1 v 1 v * M y ' S V2 

Hence = Tir^, where if = R © L 2 ((0, oo); dt) and $ : (0, oo) -> K is 

$(x) = 2" V /4 ^) © (-K/2Tr) 1/2 x 1/4 u{x)(f> t (x). (5.16) 

□ 

Remarks 5.4. (i) The condition k < in Proposition 5.3 excludes the case of the asso- 
ciated Laguerre functions w(x) = x^ l+ct ^ 2 e~ x ^ 2 L^n\x) 1 where with (n = 0, 1,2, . . .) 
are the associated Laguerre polynomials as in [14] and [7, 9.237]. 

(ii) The Laplace transforms of \fxK v {^l^fx) and y / xJ !y (2y / i : ) may be expressed in terms of 
Whittaker's functions. 

(iii) Bessel's equation may be transformed into the typical confluent hypergeometric equa- 
tion as in [8, p. 228]. 

Acknowledgement. I am grateful to the referee for pointing out reference [20]. 
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